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Abstract 

High-precision calculations of electroweak processes in light nuclei are of great im- 
portance for multiple reasons. Gerry Brown together with Dan-Olof Riska published 
in 1972 a famous calculation on radiative capture of a thermal neutron on a proton, 
and their work was a precursor of subsequent great developments in high-precision 
calculations of electroweak processes in light nuclei. The application of effective 
field theory to nuclear systems is a prominent example of these developments. We 
present an overview of it with main emphasis placed on those topics in which we 
ourselves have been involved. 



*To appear in Festschrift for Gerry Brown, ed. Sabine Lee (World Scientific, Singapore) 



1 



1 Introduction 



In 1972 Dan-Olof Riska and Gerry Brown published their celebrated calculation [1] 
on radiative capture of a thermal neutron on a proton: np — > ^d. Until their 
work appeared, the ~10 % discrepancy between the experimental and theoretical 
values of the radiative np capture rate had been a persistent puzzle. Dan-Olof 
and Gerry demonstrated that this discrepancy could be removed by including the 
contributions of the exchange currents derived by Chemtob and Rho [2], 1 and this 
was the first convincing case that pointed to the existence of the exchange currents 
in nuclei. The 1972 Riska-Brown paper was a precursor of subsequent important 
developments in high-precision calculations of electroweak processes in light nuclei. 
On the occasion of Gerry's 85th birthday, we would like to give a succinct survey of 
these developments; we shall concentrate here mostly on those topics we ourselves 
have worked on, and so this note is not meant to be a comprehensive survey of the 
field. 

High-precision calculations of electroweak processes in light nuclei are of great 
importance for multiple reasons. First of all, light nuclei, for which nuclear many- 
body-problem aspects are relatively well controlled, are expected to provide clean 
testing grounds for our understanding of the behavior of the hadronic electroweak 
currents in nuclei. In this connection, we should emphasize that the recent progress 
in formalisms and numerical techniques for solving the multi-nucleon Schrodinger 
equation has reached such a level [3] that the wave functions of low-lying levels for 
light nuclei can now be obtained practically with no approximation (once the validity 
of a given nuclear Hamiltonian is assumed). This situation frees us from the "usual" 
nuclear physics complications that result of truncating the nuclear Hilbert 

space down to certain model space (such as shell-model configurations, cluster-model 
trial functions, RPA, etc.). 

We should also mention that nuclear electroweak processes in light nuclei are 
of particular importance in connection with astrophysics and particle physics. They 
appear in many astrophysical phenomena, prominent examples being stellar evolu- 
tion (including solar nuclear fusion), big-bang nucleo-synthesis, supernova nucleo- 
synthesis, etc. Furthermore, in terrestrial experiments designed to detect astrophys- 
ical neutrinos, weak-interaction reactions on light nuclei play important roles. It 
is also to be noted that neutrino-nucleus reactions feature importantly in neutrino- 
oscillation experiments, one of the most hotly pursued frontiers in particle physics. 
Reliable information on the cross sections for these processes is obviously in great 
demand. However, for many of the processes in question, it is extremely difficult, 
or even impossible at present, to experimentally measure their cross sections. We 
must therefore totally rely on theoretical estimates, and this implies that we need 
a theoretical framework that allows us to arrive at estimates with reliable error as- 
sessments attached. Fortunately, the past two decades have witnessed significant 
progress toward this goal. 

1 see M. Rho, in this volume. 
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Recent developments in high-precision calculations of electroweak reactions in 
light nuclei encompass two lines of studies. One is the elaboration and expansion 
of the phenomenological potential model, which has come to be known as the stan- 
dard nuclear physics approach (SNPA). Another line is the successful application of 
effective field theory (EFT) to nuclear physics. In the sense that the former is phe- 
nomenological in nature whereas the latter is based on a systematic expansion that 
is formally in conformity with the tenets of EFT, these two approaches are distinct 
and, perhaps because of this, it happens rather often that practitioners of one ap- 
proach (SNPA or EFT) ignore the existence of and progress in the other approach. 
In our view, however, this does not seem an ideal or even productive situation. As 
discussed below, the forte of EFT is that it gives a clear-cut guiding principle on 
how to tame short-range physics. But, as far as long-range nuclear correlations are 
concerned, EFT does not necessarily lead to improvement over what has already 
been achieved in SNPA. And, as we will illustrate with concrete examples later, 
the high-precision calculation of an observable often requires that both short-range 
physics and long-range nuclear correlations are under control. So, at least at present, 
it appears to be a reasonable strategy to take advantage of the availability and the 
power of both SNPA and EFT while assuring, as best one can, consistency between 
the two. This viewpoint led us to propose the "more effective" effective field theory 
(MEEFT) approach, also called EFT*. It should be emphasized here that Tae-Sun 
Park has played a particularly important role in the developments of this approach. 
The nature of EFT limits us to consider only those processes which involve rela- 
tively low energies and momenta but, within this limitation, MEEFT has registered 
remarkable successes in giving reliable estimates for electroweak transition rates of 
current interest. We wish to describe here some of these results; primarily we shall 
be concerned with the following processes: 

pp-fusion: p + p — >■ d +e + +u e 

v-d reactions: u e + d — > e~+ p+p, v + d — > u+n+p 

/i-d capture: fi~ + d — > u^+n+n 

Radiatice 7r-capture: ir~ + d — y 7 + n + n. 

Tritium /3-decay: t — y 3 He + e~+ v e 

fi - 3 He capture: /i~+ 3 He — > + t 

Hep reaction 3 He +p — y 4 He +u e + e + 

Hen reaction 3 He +n — y 4 He +7 

This list may look like a haphazardly picked collection but, as will be explained 
below, they are interrelated in the context of demonstrating the important progress 
achieved in MEEFT. Specific motivations for considering these reactions will be 
described later for each case as the individual processes are discussed. Just to set 
the stage, however, we make some remarks here on the following three processes: 
pp-fusion, Hep and the vd reaction. 
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The pp fusion reaction, pp — > de + u e , is the primary solar nuclear reaction 
that essentially controls the burning rate of the sun, and hence the exact value of 
its cross section is a crucial input for any further developments of solar models. 
According to the latest review [4], this rate is required to be known with better 
than ~1 % precision in order to give impact to the existing models. Meanwhile, the 
Hep reaction, 3 He+p — > A He + e + + v e is important in a different context. Although 
Hep plays only a very minor role in the pp-chain solar burning, it generates highest- 
energy solar neutrinos whose spectrum extends beyond the maximum energy of 
the 8 B neutrinos. So, even though the flux of the Hep neutrinos is small, it can 
distort the higher end of the 8 B neutrino spectrum [5], and this distortion can affect 
the interpretation of the results of a Super-Kamiokande experiment [6]. On the 
theoretical side, as explained in more detail below, Hep represents a unique case in 
which the leading-order contribution is drastically suppressed. This implies that Hep 
offers a good place for investigating higher order terms in EFT expansion. As for 
the v-d reactions, their primary importance is connected to the Sudbury Neutrino 
Observatory (SNO) experiments, in which the solar neutrino flux was monitored 
using a heavy water Cerenkov counter. As is well known, the SNO experiments [7] 
have established that the total solar neutrino flux (summed over all neutrino flavors) 
agrees with the prediction of the standard solar model [8], whereas the electron 
neutrino flux from the sun is significantly smaller than the total solar neutrino flux. 
These results of the SNO experiments have given "smoking-gun" evidence for the 
transmutation of solar electron neutrinos into neutrinos of other flavors. Obviously, 
a precise knowledge of the u-d reaction cross sections is of primary importance in 
interpreting the SNO data. 

Before going into the description of achievements obtained in nuclear EFT, 
it seems useful to take a quick look at the latest lattice QCD calculation [9] of the 
nucleon weak form factors defined as 



<n\V tl \p> = u r 
<n\A fl \p> = u r 



2m N 



Up , 



u p , (1) 



where q = p n — p p . The calculation [9] was carried out for the u/d quark mass 
that corresponded to m n = 0.33 GeV ~ 0.67 GeV, and the form factors were cal- 
culated for 0.2 GeV 2 < q 2 < 0.75 GeV 2 . Here is a partial summary of the results: 
(1) Fy(q 2 ) is described by the conventional dipole form, but the dipole mass Mi is 
found to be 27 % smaller than the empirical value; (2) FM{q 2 ) is described by the 
dipole form, but the corresponding dipole mass M2 is found to be 26 % smaller than 
the empirical value; (3) gA/gv — 119, which is 7 % smaller than the experimen- 
tal value (g A /gv)e xp = 1-2695 ± 0.0029 (PDG08); (4) F A (q 2 ) is found to have the 
dipole form but the axial dipole mass Ma is found to be 37 % larger than M^ xp ; 
(^) ^2m^i^(g2) ' = &(q 2 ) is found to be constant in conformity with pion-pole dom- 
inance, but a = 0.8 ~ 0.9 (instead of a — 1 expected from PCAC). These features 
lead us to expect that it will be a while before nuclear elecroweak processes can be 
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calculated with required precision by lattice QCD. 



2 Theoretical frameworks 

• Standard nuclear physics approach (SNPA) 

The phenomenological potential picture has been highly successful in describ- 
ing many different nuclear phenomena. In this picture an A-nucleon system is 
described by a Hamiltonian of the form 

A A A 

H = E^ + E^ + E Vijk + ---, (2) 

i i<j i<j<k 

where t{ is the kinetic energy of the i-th nucleon, is a phenomenological two-body 
potential, Vijk is a phenomenological three-body potential, and so on. Once H is 
specified, the A-body nuclear wave function |\&> is obtained by solving exactly or 
essentially exactly the Schrodinger equation 

H\^>= E\^> . (3) 

There is much freedom in choosing possible forms of apart from a well-established 
requirement that, as the inter-nucleon distance becomes sufficiently large, 
should approach the one-pion exchange Yukawa potential. For the model-dependent 
short-range part of Vij, one assumes certain functional forms and fix the parame- 
ters appearing therein by demanding that the solutions of eq.(3) for the A=2 case 
reproduce the nucleon-nucleon scattering data (typically up to the pion-production 
threshold energy) as well as the deuteron properties. There exist several so-called 
high-precision phenomenological N-N potentials that can reproduce all the existing 
two- nucleon data with normalized x 2 values close to 1. These potentials differ sig- 
nificantly in the ways they parameterize short-range physics and, as a consequence, 
they exhibit substantial difference in their off-shell behavior. An important ques- 
tion therefore is to what extent the model-dependent short-range part affects the 
reliabilities of theoretical values of various observables in A=2 and higher systems 
calculated with these phenomenological potentials. 

Nuclear responses to external electroweak probes in normal circumstances are 
given, to good approximation, by one-body terms, which are also called the impulse 
approximation (IA) terms. To achieve higher accuracy, however, we must also con- 
sider exchange current (EXC) terms, which represent the contributions of nuclear 
responses involving two or more nucleons. If for some reason the IA contributions 
are suppressed, it becomes particularly important to consider the EXC contribu- 
tions. These exchange currents are usually taken to be two-body operators, and 
they are derived from one-boson exchange diagrams; the vertices featuring in the 
relevant diagrams are determined to satisfy the requirements of low-energy theorems 
and current algebra [2, 10]. A formalism based on this picture is commonly called 
the standard nuclear physics approach (SNPA). (This is also called a potential model 
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in the literature.) Schematically, an electroweak nuclear matrix element in SNPA is 
given by 

M s f ? PA =<*f\T,Ot+T,°tm+ Ee>wl*i>, (4) 

l £<m £<m<n 

where the initial (final) nuclear wave function, (^/), is a solution of eq.(3), and 
C?£, Oe m and Oi mn are, respectively, one-body, two-body and three-body transition 
operators for a given electroweak process. SNPA has been used extensively for 
describing electroweak processes in light nuclei, and general good agreement found 
between theory and experiment [3] strongly suggests that SNPA captures much of 
the physics involved. 

• Effective field theory (EFT) 

Even though SNPA has been scoring great successes in describing a vast va- 
riety of data, it is still important from the fundamental point of view to raise the 
following issues. First, given the fact that hadrons and hadronic systems (such as 
nuclei) are governed by quantum chromodynamics (QCD), it is desirable to estab- 
lish connections between SNPA with QCD. In this connection we note that, whereas 
chiral symmetry is known to be a fundamental symmetry of QCD, the formulation 
of SNPA is largely disjoint from this symmetry. Secondly, in introducing a phe- 
nomenological SNPA Lagrangian for describing the nuclear interaction and nuclear 
responses to the electroweak currents, it is not clear whether there is any identifi- 
able expansion parameter that helps us to control the possible forms of terms in the 
Lagrangian and that gives a general measure of errors in our calculation. 

An EFT-based approach [11] is a powerful tool to address these issues, and it 
has been studied with great intensity; for reviews, see [12]- [16]. The basic idea of 
EFT is simple. To describe phenomena characterized by a typical energy- momentum 
scale Q, we expect that we can exclude from our Lagrangian those degrees of freedom 
that belong to energy-momentum scales much higher than Q. So we introduce a 
cut-off scale A that is sufficiently larger than Q and we categorize our fields (to be 
generically represented by 0) into two groups: high-frequency fields H and low- 
frequency fields L . By integrating out H , we arrive at an effective Lagrangian that 
only contains L as explicit dynamical variables. The effective Lagrangian £ eff is 
related to the original Lagrangian C as 

j[d(f ) }e [ I d4xC ^ = j[d(f )H }[d ( j) L }e i I d4xC ^' M = j[dcf)^e^ dixCM . (5) 

£ e ff defined by eq.(5) inherits the symmetries (and the pattern of symmetry breaking, 
if there is any) of the original Lagrangian C. It can be shown that £ c fr is given by the 
sum of all possible monomials of L and their derivatives that are consistent with the 
symmetry requirements dictated by C Because a term involving n derivatives scales 
like (Q/A) n , the terms in £ eff can be organized into a perturbative series wherein 
Q/A serves as an expansion parameter. The coefficients of terms in this expansion 
are called low-energy constants (LECs). If all the LEC's up to a specified order n 
can be fixed either from theory or from fitting to the experimental values of relevant 
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observables, C e s serves as a complete (and hence model-independent) Lagrangian 
to the given order of expansion. The results obtained are ensured to have accuracy 
of order (Q/A) n+1 . In applying EFT to nuclear systems, the underlying Lagrangian 
C is the QCD Lagrangian Cqcd, whereas, for the typical nuclear physics energy- 
momentum scale Q <C A x ~ 1 GeV, the effective degrees of freedom that should 
feature in £ cff are hadrons rather than the quarks and gluons. It is not trivial 
to apply the formal definition in eq.(5) to derive C c g written in terms of hadrons 
starting from Cqcd] this is because the hadrons cannot be simply identified with 
the low-frequency field, <f>L in eq.(5). At present, the best one could do is to resort to 
symmetry considerations and the above-mentioned expansion scheme in the spirit 
of Weinberg's folk theorem [17]. Here chiral symmetry plays an important role. 
We know that chiral symmetry is spontaneously broken, generating the pions as 
Nambu-Goldstone bosons. 2 This feature can be taken into account by assigning 
suitable chiral transformation properties to the Goldstone bosons and writing down 
all possible chiral-invariant terms up to a specified chiral order [18]. The above 
consideration presupposes exact chiral symmetry in Cqcd- In reality, however, Cqcd 
contains small but finite quark mass terms, which explicitly break chiral symmetry 
and lead to a non- vanishing value of the pion mass m n . Again, there is a well- 
defined method to determine what terms are needed in the Goldstone boson sector 
to represent the effect of explicit chiral symmetry breaking [18]. These considerations 
lead to an EFT called chiral perturbation theory (%PT) [17, 19]. The successes of 
xPT in the meson sector are well documented; see, e.g., [12]. 

A serious problem we encounter in extending to the nucleon sector is 
that, as the nucleon mass m N is comparable to the cut-off scale A x , a simple ap- 
plication of expansion in Q/A does not work. This difficulty can be circumvented 
by employing heavy-baryon chiral perturbation theory (HBxPT), which essentially 
consists in shifting the reference point of the nucleon energy from to m N and in 
integrating out the small component of the nucleon field as well as the anti-nucleonic 
degrees of freedom. Thus an effective Lagrangian in HB^PT contains as explicit de- 
grees of freedom the pions and the large components of the redefined nucleon field. 
HB^PT has as expansion parameters Q/A x , m 7T /A x and Q/m N . As m N A x , it is 
convenient to combine chiral and heavy-baryon expansions and introduce the chiral 
index v defined by v = d + (n/2) — 2. Here n is the number of fermion lines that 
participate in a given vertex, and d is the number of derivatives (with m n counted 
as one derivative). A similar power counting scheme can be introduced for Feynman 
diagrams as well [11]. The contribution of a Feynman diagram that contains Na 
nucleons, Ne external fields, L loops and Nc disjoint parts is shown to scale like 
(Q/A) u , where the chiral index v is defined as 

u = 2L + 2(N c -l) + 2-(N A + N E ) + J2^i, (6) 

i 

with the summation running over all the vertices contained in the Feynman diagram. 
HBxPT has been applied with great success to the one-nucleon sector [12]. 

2 We limit ourselves here to SU(2)xSU(2) chiral symmetry. 
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However, HB%PT cannot be applied in a straightforward manner to nuclei, 
which contain more than one nucleon. The reason is that nuclei involve very low- 
lying excited states, and the presence of this small energy scale spoils the original 
counting rule [11]. A method to avoid this difficulty was proposed by Weinberg [11]. 
Classify Feynman diagrams into two categories, irreducible and reducible diagrams. 
Irreducible diagrams are those in which every intermediate state has at least one 
meson in flight; all others are classified as reducible diagrams. One then applies 
the above chiral counting rules only to irreducible diagrams. The contributions of 
all the two-body irreducible diagrams (up to a specified chiral order) are treated as 
an effective potential (to be denoted by V FFT ) acting on nuclear wave functions. 
Meanwhile, the contributions of reducible diagrams can be incorporated [11] by 
solving the Schrodinger equation 

H EFT \m> EFT = E\V> EFT , (7) 

where 

H EFT = X> + tv* FT + E V*r + --- , (8) 

i i<j i<j<k 

We may refer to this two-step procedure as nuclear %PT or, more specifically, nuclear 
%PT in the Weinberg scheme. 3 

To apply nuclear xPT to a nuclear electroweak process, we derive a nuclear 
transition operator T from a complete set of all the irreducible diagrams (up to a 
given chiral order v) involving the relevant external current (s) [20, 21, 22]. Thus 

A A 

r FT = J2 °T T + E °T + ■■■ ( 9 ) 

i i<j 

To preserve consistency in chiral counting, the matrix element of 7~ EFT must be 
calculated with the use of nuclear wave functions generated with the use of nuclear 
interactions that represent all the irreducible A-nucleon diagrams up to v-th order. 
Thus a transition matrix in nuclear EFT is given by 

A A 

m^j t = < ^ FFT it 551 ^ |\i/f FT >=< ^ FFT | ®T T + E ^y FT + " ' l^f FT> > ( 10 ) 

i i<j 

where the superscript, "EFT", means that the relevant quantities are calculated 
according to EFT as described above. If this program is carried out rigorously, it 
would constitute an ab initio calculation. Let us call it "AB-EFT" . It is to be noted 
that in EFT we know exactly at what chiral order three-body operators start to 
contribute to T, and that, to chiral orders relevant to most applications described 
below, there is no need for three-body operators. With this understanding, we have 
kept only one- and two-body operators in eq.(10). The unambiguous classification 
of transition operators according to their chiral orders is a great advantage of EFT, 
which is missing in eq.(4). 



3 This is often called the A-counting scheme [14]. 
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We should point out that there exists an alternative form of nuclear EFT called 
the PDS scheme, proposed by Kaplan, Savage and Wise [23], which uses a counting 
scheme (often called Q-counting) different from the Weinberg scheme. Although 
many important results have been obtained in the PDS scheme (for a review, see 
e.g. [15]), we concentrate here primarily on the Weinberg scheme. 

We also remark that, for very low energy nuclear phenomena whose typical 
energy-momentum scale is Q <C m n , even the pions can be treated as "heavy" 
particles and can be eliminated from C e s, with the choice of A m n . This leads to 
pionless EFT [23, 24, 25]. Here there is no manifest chiral symmetry since there is 
no explicit pion. In this regime the NN interactions and electroweak currents are 
described by point-like contact terms. 

• Hybrid EFT 

In the above we emphasized the formal merits of nuclear EFT. In actual calcu- 
lations, however, we need to consider the following two points. First, it is still a big 
challenge to generate, strictly within the EFT framework, nuclear wave functions 
the accuracy of which is comparable to that of SNPA wave functions. Secondly, 
whereas a chiral Lagrangian, £ eff , is definite only when the values of all the rele- 
vant LECs are fixed, this condition cannot be readily satisfied in many cases. A 
pragmatic solution to the first difficulty [26]- [31] is to use in eq.(10) wave functions 
obtained in SNPA; this eclectic approach is called hybrid EFT. A nuclear transition 
matrix element in hybrid EFT is given by 

M)r EFT = < ^ npa i £ of FT + £ oz t \^ npa > , (ii) 

Because, as mentioned, the NN interactions that generate SNPA wave functions 
reproduce accurately the two-nucleon data, the adoption of eq.(ll) is almost equiv- 
alent to using the empirical data themselves to control the initial and final nuclear 
wave functions. From a purely formal point of view, hybrid EFT may be deemed 
as a "regression". But, if our goal is to obtain a transition matrix element as accu- 
rately as possible with the maximum help of available empirical input, hybrid EFT 
does have a legitimate status. It offers more flexibility and predictive power than 
the AB-EFT, and it can be highly accurate so long as the off-shell problem and the 
contributions of three-body (and higher-body) terms are properly addressed. The 
issue of possible unknown LECs will be discussed in connection with MEEFT in the 
next subsection. 

• MEEFT 

A great merit of hybrid EFT is that it can be used for complex nuclei (A = 3, 
4, ...) with essentially the same accuracy and ease as for the A=2 system. 4 We 

4 Hcre we are ignoring "purely technical" complications that can grow in actual numerical cal- 
culations for higher- A systems. 
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reemphasize in this connection that the contributions of transition operators involv- 
ing three or more nucleons, which can in principle contribute in A-nucleon systems 
(A > 3), are intrinsically suppressed according to chiral counting, and hence, up 
to a certain chiral order, a transition operator in an A-nucleon system consists of 
the same EFT-based 1-body and 2-body terms as used for the two-nucleon system. 
Then, since SNPA provides high-quality wave functions for the A-nucleon system, 
one can calculate J\A t Jf b ~ EFT with precision comparable to that for the corresponding 
two-nucleon case. 

In most cases the one-body operator, Of FT , is free from unknown LECs. So 
we may concentrate on the two-body operator, Of^ T . Suppose that OfJ^ T under 
consideration contains an LEC (call it k) that cannot be determined with the use of 
A=2 data alone. It can happen that an observable (call it f2) in a A-body system 
(A > 3) is sensitive to k and that the experimental value of Q is known with sufficient 
accuracy. We then can determine k by calculating Ai f Jf b ~ EFT responsible for Q and 
adjusting k to reproduce the experimental value of Q. Once k is determined this 
way, we can make predictions for any other observables for any other nuclear systems 
that are controlled by the same transition operators. When hybrid EFT is used in 
this manner, it is called more effective effective field theory (MEEFT) [30, 31]. 5 

MEEFT is a powerful formalism for correlating various observables in different 
nuclei, using the transition operators controlled by EFT. A further practical advan- 
tage of MEEFT is that, since correlating the observables in neighboring nuclei is 
expected to serve as an additional renormalization, the possible effects of higher chi- 
ral order terms and/or off-shell ambiguities can be significantly suppressed with the 
use of MEEFT. We will come back to this later, when we discuss concrete examples. 



3 Examples 

• pp-fusion, vol reactions and fid capture 

A common feature of these reactions is that a precise knowledge of the Gamow- 
Teller (GT) transition matrix elements is required in estimating their cross sections. 
We show here, following Refs. [31], that the idea of MEEFT can be used very neatly 
for this group of reactions. The 1-body IA operators for the GT transition can be 
fixed unambiguously from the available 1-body data. As for the 2-body operators, 
to next-to-next-to-next-to-leading order (N 3 LO) in chiral counting, there appears 
one unknown LEC that at present cannot be determined from data for the A=2 
systems. This unknown LEC, denoted by d R in [27], parametrizes the strength of 
contact-type four-nucleon coupling to the axial current. Park et al. [31] noticed that 
the same LEC, d^, also features as the only unknown parameter in the calculation of 
the tritium /3-decay rate F^, and they proposed to use MEEFT to place a constraint 

5 It is also called EFT*. MEEFT should be distinguished from an earlier naive hybrid EFT 
model wherein the short-range terms were dropped altogether using an intuitive argument based 
on short-range NN repulsion. 
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on dji from the experimental information on Fg. Since the empirical value of rjg 
is known with high precision, and since the accurate wave functions of 3 H and 3 He 
are available from SNPA [32], we can determine da with sufficient accuracy for our 
purposes. Once the value of d,R is determined this way, we can carry out parameter- 
free MEEFT calculations for pp-fusion [31], v-d reactions [33], and /i-d capture. (An 
on-going attempt to deduce d,R from /i-d capture will be discussed later.) 

We should mention the important role of momentum cutoff in MEEFT. As 
emphasized before, the effective Lagrangian £ e // is, by construction, valid only be- 
low the specified cutoff scale A. Needless to say, this basic constraint should be 
respected in our nuclear EFT calculations, and for that we must make sure that nu- 
clear intermediate states involved in the computation of eq.(10) do not get out of this 
constrained world. Since, however, the use of \/ phcnom introduces high momentum 
components above Aq CD , we need to introduce momentum cutoff at A NN to elim- 
inate high-momentum components. It is reasonable to implement this constraint 
by requiring that the two-nucleon relative momentum should be smaller than Ann- 
Park et al. [31] used a Gaussian cutoff function proportional to exp(— p 1 / A 2 NN ) but 
its detailed form should not be too important. As a reasonable range of the value 
of A NN we may choose: 500 MeV < A^tv^ 800 MeV, where the lower bound is 
dictated by the requirement that Ann should be sufficiently larger than the pion 
mass (to fully accommodate pion physics), while the upper bound reflects the fact 
that our EFT is devoid of the p meson. Mismatch between v|/ EFT and \[/ SNPA only 
affects short-distance behavior, and we expect that A^Ar-independence is a measure 
of model independence of an MEEFT calculation. 

For a given value of A NN within the above range, d,R is tuned to reproduce 
and then the rates for pp-fusion, the vd reactions and fid capture are calcu- 
lated. Before giving an account of the individual results, we should point out a 
notable common feature. Although the optimal value of <1r varies significantly as 
a function of A N n, the observables (in our case the above three reaction rates) ex- 
hibit remarkable stability against the variation of A NN (within the above-discussed 
physically reasonable range). This stability may be taken as an indication that the 
use of MEEFT for inter-correlating the observables in neighboring nuclei effectively 
renormalizes various effects, such as the contributions of higher-chiral order terms, 
mismatch between the SNPA and EFT wave functions, etc. This stability is essential 
in order for MEEFT to maintain its predictive power. 



3.1 The ^-fusion S-factor 

The .S-factor is related to the cross section a as 

a(E) = ^-exp{-2n V (E)}, 

where E is the incident energy in the center-of-mass system, and rj(E) = Zx ^ff is 
the Coulomb penetration factor. The .S-factor for pp-fusion at threshold is related 
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to the transition matrix element M as [4] 

S pp (0) = 6n 2 m p a 1d2 \(g A /gv) 2 , ,J™ \M\ 2 (12) 

where 7 = (2fxB ( i) 1 ^ 2 = 0.2316 fm -1 is the deuteron binding wave number, and 
is the phase space factor including radiative corrections. As mentioned, the solar 
model and stellar evolution theory require 1% precision in S pp [4]. 

It is useful to first recall what was achieved in SNPA. A benchmark calculation 
of S pp in SNPA was carried out by Schiavilla et al. [34]. These authors took full 
advantage of the fact, first pointed out in Ref. [35], that the Gamow- Teller (GT) 
matrix elements for pp-fusion and tritium /3-decay are closely correlated. For each 
of the two-body GT transition operators with different r-dependences (r = relative 
distance of the two nucleons), the transition density for pp-fusion is found to be 
accurately proportional to that for tritium /3-decay. This means that a single multi- 
plicative constant can match them all. Schiavilla et al. [34] fine-tuned the iVA-axial 
coupling constant, which featured in a dominant two-body axial current diagram 
(A-particle excitation diagram), and which was only loosely controlled from other 
sources of information. After this fine-tuning, Spp was calculated with estimated 
uncertainties of ~0.3 %; here the error estimate came from the range of variation in 
S pp (0) for the five different high-precision phenomenological potentials used in [34]. 
A further examination of this error estimation will be described below. 

Park et al. [31] carried out an MEEFT calculation of S pp (0) and obtained 

Spp(0) = 3.94 x (1 ± 0.008) x 10" 25 MeV b . (13) 

S P p(0) is found to change only by ~0.1% against changes in Ann, ensuring the 
robustness of the MEEFT calculation in this respect. The convergence pattern 
in chiral expansion for the pp-fusion transition amplitude is reported to be [31]: 
M ~ (1 + + 0.1% + 0.9% + 0.4%). It is noted that the 2nd and 3rd order terms 
are accidentally small, but that the 4th and 5th order are roughly in conformity 
with expansion in m 7T /A x |. The 3-body currents contribute at the sixth order. 
Uncertainties due to the higher order contributions estimated from the size of the 
5th-order contribution are 0.4 % (in amplitude), or 0.8 % (in probability). 6 The 
MEEFT result, eq.(13), is consistent with that obtained in SNPA [34]. Meanwhile, 
the fact that MEEFT allows us to make error estimation is a notable advantage over 
SNPA. It is worth emphasizing that the accuracy achieved in eq.(13) represents an 
improvement by a factor of ~5 over the previous results based on a naive hybrid 
EFT [27]. 

The calculation of S pp in pionless EFT was performed up to NNLO by Kong 
and Ravndal [36], and up to fifth order by Butler and Chen [37]; see also Ando 

6 In [31], the relative error of 0.4 % was assigned to S pp (0). However, it has been found subse- 
quently that the 0.4 % uncertainty should have been assigned to the transition amplitude, not to 
the transition probability; we therefore double the error estimate quoted in [31] and adopt 0.008 
in eq.(13). 
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et al. [38]. Here again, up to the order considered, there is only one unknown 
LEC, denoted by Li ^, which (analogously to d R ) parameterizes the strength of 
four-nucleon contact axial-vector coupling. The results of Refs. [36, 37] are: (i) 
the central value of 5^(0) is consistent with those obtained in SNPA and MEEFT; 
(ii) the theoretical uncertainty in S pp (0) is ~3 %, which is significantly larger than 
the results obtained in SNPA and MEEFT. One thing to be noted here is that the 
PDS scheme [23] (used in pionless EFT) adopts an expansion scheme for transition 
amplitudes themselves, without employing the concepts of the potential or wave 
functions. In certain contexts this feature may be an advantage, but its disadvantage 
in the present context is that one cannot readily relate the transition matrix elements 
for an A-nucleon system with those for the neighboring nuclei; in PDS, each nuclear 
system requires a separate parameterization. Therefore, L\a in pionless EFT needs 
to be determined from the empirical value of a 2-body observable, which however is 
at present available only with large experimental error. This situation is reflected 
in the rather large 3 % uncertainty in S pp {$$) calculated in pionless EFT. 

In a very recent review [4], there has been an attempt to deduce the most 
updated estimate for S pp (0), and the recommended value is 

Spp(0) = 4.01 x (1 ± 0.009) x 10~ 24 MeVb. (14) 

This estimate arises as follows. In eq.(12), the matrix element M is assigned the 
value corresponding to the central value obtained in SNPA [34] and MEEFT [31]; 
these two approaches give exactly the same value. For the other inputs that enter 
into eq.(12), the most updated values are used; some of them slightly differ from 
the values used in the earlier publications. The relative uncertainty of 0.009 in 
eq.(14) is dominated by the 0.8 % uncertainty in \M\ 2 estimated from the MEEFT 
result [31] and the 0.5 % uncertainty coming from the experimental error in the 
nn scattering length a nn . These uncertainties and other smaller uncertainties are 
added in quadrature. Eq.(14) indicates that S pp (0) now can be estimated within 1% 
precision. 

3.2 v-d reactions 

As mentioned, v-d reactions are important in connection with the SNO experiments 
on neutrino oscillations. Within SNPA a detailed calculation of the v-d cross sec- 
tions, a(vd), was carried out in Refs. [39, 40, 41]. Nakamura et al. [41], following 
Schiavilla et a/.'s approach [34], used the two-body current whose strength was 
fine-tuned to reproduce the tritium /3-decay rate. With the use of this fine-tuned 
two-body exchange current, Nakamura et al. calculated a{vd) for four types of re- 
actions: charged-current and neutral-current reactions on d for both neutrino and 
anti-neutrino cases. Meanwhile, Butler, Chen and Kong (BCK) [42] carried out a 
pionless-EFT calculation of the v-d cross sections. The results obtained agree with 
those of Nakamura et al. [41] in the following sense. A pionless-EFT calculation, 
as mentioned in the discussion of pp-fusion, involves one unknown LEC, L^, which 
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represents the strength of a four-nucleon axial-current coupling term. BCK deter- 
mined by requiring that the vd cross sections of NETAL be reproduced by their 
pionless-EFT calculation. With the value of L\ A fine-tuned this way, the <r(z/<i)'s 
obtained by BCK show a perfect agreement with those of Nakamura et al. for all 
the four reactions mentioned above, and for the entire solar neutrino energy range, 
E„< 20 MeV. Moreover, the optimal value, L 1A = 5.6 fm 3 , found by BCK [42] is 
consistent with the order of magnitude of Lia expected from the naturalness ar- 
gument (based on a dimensional analysis), |Lia| < 6fm 3 . The fact that an EFT 
calculation (with one parameter fine-tuned) reproduces the results of SNPA very 
well strongly suggests the robustness of the SNPA results for a(vd). 

Even though it is reassuring that the v-d cross sections calculated in SNPA 
and pionless-EFT agree with each other (in the sense explained above), it is desirable 
to carry out an EFT calculation that is free from any adjustable LEC. Fortunately, 
MEEFT allows us to carry out an EFT-controlled parameter-free calculation of the 
v-d cross sections, because the only LEC appearing therein, du, has already been 
determined from [31]. Ando et al. [33] took advantage of this situation and carried 
out an MEEFT calculation of a(vd). The results were found to agree within 1% 
with a(vdys obtained in SNPA [41]. These results show that the v-d cross sections 
used in interpreting the SNO experiments [7] are reliable at the 1% precision level. 

There have been subsequent attempts to calibrate L\a from reactor z/-deuteron 
breakup reactions [43], and also from a "self-calibrating" analysis of the SNO data [44] 
The a(z/<i)'s corresponding to the range of L±a obtained in these analyses are con- 
sistent with the a(z/<i)'s obtained in SNPA and MEEFT. However, the resulting 
constraints on L\a are not stringent. We have mentioned that both L\a and d,R 
represent the strength of axial-current-four-nucleon contact coupling. It is to be 
noted, however, that L\a belongs to pion-less EFT, while dn to pion-ful EFT. In 
the pion-ful EFT, because of the strong tensor force, the exchange current involving 
the deuteron d-state is important, and the s-wave exchange current arising from the 
d R term is separate from this tensor- force effect. By contrast, in the pion-less EFT, 
the explicit rf-wave term is a higher-order correction, and hence the s-wave L\a 
term must subsume the strong tensor-force contributions. It would be informative 
to investigate the relation between L\a and c/r from this perspective. 

• Off-shell effects 

In introducing hybrid EFT, we have replaced |\P > EFT for the initial and final 
nuclear states in eq.(10) with the corresponding |\P > 5JVPA ; see eq.(ll). This replace- 
ment may bring in a certain degree of model dependence, called the off-shell effect, 
because the phenomenological NN interactions are constrained only by the on-shell 
two-nucleon observables. 7 This off-shell effect, however, is expected to be small for 
the reactions under consideration, since they involve low momentum transfers and 

7 In a consistent theory, physical observables are independent of field transformations that lead to 
different off-shell behaviors, and therefore the so-called off-shell effect is not really a physical effect. 
In an approximate theory, observables may exhibit superficial dependence on off-shell behavior, 
and it is customary to refer to this dependence as an off-shell effect. 
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hence are not extremely sensitive to the short-range behavior of the nuclear wave 
functions. One way to quantify this expectation is to compare a two-nucleon rela- 
tive wave function generated by the phenomenological potential with that generated 
by an EFT-motivated potential. Phillips and Cohen [45] made such a comparison 
in their analysis of the 1-body operators responsible for electron-deuteron Comp- 
ton scattering, and showed that a hybrid EFT should work well up to momentum 
transfer 700 MeV. A similar conclusion is expected to hold for a two-body operator, 
so long as its radial behavior is duly "smeared-out" reflecting a finite momentum 
cutoff. Thus, hybrid EFT as applied to low energy phenomena is expected to be 
practically free from the off-shell ambiguities. The off-shell effect should be even less 
significant in MEEFT, wherein an additional "effective" renormalization is likely to 
be at work (see subsection 3.1). 

We now discuss briefly another interesting development, due to Tom Kuo and 
his colleagues [46], which can shed much light on the reliability of a hybrid EFT 
or MEEFT calculation. As mentioned, a "realistic phenomenological" nuclear in- 
teraction, Vij in eq.(2), is determined by solving the Schrodinger equation, eq.(3), 
for the A=2 system and fitting the results to the full set of two-nucleon data up 
to the pion production threshold energy. So, physically, should reside in a mo- 
mentum regime below a certain cutoff, A c . In the conventional treatment, however, 
the existence of this cutoff scale is ignored, and eq.(3) is solved, allowing the en- 
tire momentum range to participate. Kuo et al. proposed to construct an effective 
low-momentum potential Vi ow -k by eliminating (or integrating out) from the 
momentum components higher than A c , and calculated V/ otl) _fc's corresponding to 
many well-established examples of V^'s. Remarkably, it was found that all these 
Vi ow -kS give identical results for the half-off-shell T-matrices, even though the ways 
short-range physics is encoded in these V^'s are highly diverse. This implies that 
the Vimu-kS are free from the off-shell ambiguities, and therefore the use of V/ O10 __fc's 
is essentially equivalent to employing V FFT (that appeared in eq.(8)), which by 
construction should be model- independent. Now, as mentioned, our MEEFT calcu- 
lation has a momentum-cutoff regulator, and this essentially ensures that the matrix 
element, Ai l ^ b ~ EFT ', in eq.(ll) is only sensitive to the half-off-shell T-matrices that 
are controlled by Viaw-k instead of V^. Therefore, we can expect that the MEEFT 
results reported here are essentially free from the off-shell ambiguities. 

Recently, Mosconi et al. [47] have compared vd cross sections calculated for 
various models that differ in the way the "nuclear PCAC requirement" is imple- 
mented, and concluded from this comparison that the cr(i/rf)'s obtained in SNPA, 
MEEFT and pionless EFT have uncertainties as large as 2-3 %. Mosconi et al. 
have reached this conclusion by examining the scatter of unconstrained calculations 
of a(isd). However, all state-of-the-art calculations use to reduce two-body cur- 
rent and other uncertainties. Once this constraint is imposed, the scatter in the 
calculated values of a{vd) should be significantly reduced. 
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3.3 fid capture and MuSun experiment 

Although it is likely that the determination of dp from is good enough for practical 
purposes, it is worthwhile to study a possibility to fix dp with the use of an observable 
belonging to the A=2 systems. A promising candidate is the /^-capture process, 
fi~ + d — > u^ + n + n. Even though rather large energy-momentum transfers caused 
by the disappearance of a fi~ seem to make the applicability of EFT here a delicate 
issue, we can show [48] that, as far as the hadron sector is concerned, fi-d capture is 
in fact a reasonably "gentle" process. This is because: (1) the z/ M carries away most 
of the energy, and (2) there is a large enhancement of the transition amplitude in a 
kinematic region where the relative motion of the final two nucleons is low enough to 
justify the use of EFT. Thus fid capture can be used for controlling d R , if the quality 
of experimental data improves sufficiently The MuSun experiment at PSI (proposed 
in 2008) [49] aims to determine the fid capture rate Y^d with ~1.5 % accuracy; for 
a review, see [50]. This will enable us to determine d 11 within A=2 system. In this 
connection, it is important to point out that gp = Fa(o 2 = — 0.88m^) is totally 
under control by now. Previously, the main goal of /i-nucleus capture experiments 
was to obtain information on gp given the fact that fi capture goes much faster on 
nuclei (with high Z) than on a proton. This goal, however, has been perennially 
elusive because of nuclear physics complications. It is therefore highly noteworthy 
that a HB^PT calculation [51] gives a robust prediction, g P = 8.26 ± 0.23, and 
that this value is in good agreement with the experimental value, gp = 7.7 ± 1.1, 
obtained in the MuCap Experiment at PSI [52]. The fact that the 1-body current 
is well known allows us to concentrate on the 2-body part. 

Recent SNPA calculations of fi-d capture can be found in Refs. [53, 54, 55]; 
According to Tatara et al. [53], the hyperfine-doublet capture rate Yd (by far the 
dominant capture channel) is Yd = 397 — 400 s -1 , and the results are stable against 
the use of different NN potentials available at the time of their work. The lat- 
est calculation by Ricci et al. [55] reports Y d = 416 s _1 , and that there is non- 
negligible iViV-potential dependence between the earlier and new-generation real- 
istic phenomenological potentials. Thus the results obtained within SNPA exhibit 
some scatter, and although Ricci et al. offer rather detailed discussion of its pos- 
sible origin, further investigation seems needed. An MEEFT calculation of Y d was 
carried out by Ando et al. [48], yielding the value Yd = 386 s -1 with a high degree 
of stability against variations in A^n- A pionless EFT calculation was performed 
by Chen et al. [56] but this approach is reliable only for limited kinematical regions, 
because of the very low momentum cutoff inherent with pionless EFT. 

Comparison with experiment is complicated by the existence of two barely 
overlapping experimental results: T^ xp = 470 ± 29 s -1 [57], and T^ xp = 409 ± 
40 s -1 [58]. All the above-mentioned theoretical values are consistent with the value 
of T^ xp reported in [58], and disagree with that given in [57]. 

Since the MuSun experiment envisages to measure Yd with 1.5 % precision, it 
is important to match the accuracy of theoretical prediction. Most recently, a highly 
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elaborate calculation of T d has been performed in both SNPA and MEEFT by Mar- 
cucci et al. [59] 8 . The nuclear wave functions are generated from the Argonne Vis 
potential (for SNPA) or from a chiral N3LO two-nucleon potential [60] (for MEEFT). 
As in the pp-fusion calculation, one parameter needs to be fixed (the iV-to-A axial 
coupling constant in SNPA or cLr in MEEFT). This is fixed by reproducing the tri- 
tium /3-decay rate. The model dependence with respect to the adopted interactions 
and current (including cutoff dependence in the case of MEEFT) is found to be 
weak. Marcucci et al. give 

Y d = 392 ± 2.3 s" 1 . (15) 

It is interesting to note that this result is close to that obtained by Ando et al. [48] . 
We also remark that what Marcucci et al. refer to as MEEFT is actually "less 
hybrid" in nature than the MEEFT used by Ando et al. in that the former employs 
the chiral N3LO two-nucleon potential instead of a phenomenological NN potential. 
When the MuSun result becomes available, one can use it in two ways - compare 
it with the theoretical prediction in eq.(15) based on T^, or use it to determine dft 
from the A=2 system (without referring to the A=3 system). 

3.4 Hep and Hen reactions 
• Hep reaction 

As mentioned, the Hep reaction, 3 He + p — > 4 He + v e + e + , produces solar 
neutrinos of the highest maximum energy (E™ ax = 18.8 MeV), although the flux is 
extremely low. An accurate estimation of this cross section is a particularly chal- 
lenging task because: (1) the contribution of the leading-order 1-body GT operator 
is highly suppressed due to the approximate wave function orthogonality [61], and 
(2) there is a strong cancellation between the 1-body and 2-body GT matrix el- 
ements [35, 32]. Reflecting this situation, estimates of 5 , H ep (0) changed by orders 
of magnitude from Salpeter's original value, SHep(0)s a ipeter = 630 x 10~ 20 keV-b to 
modern values; for a review, see [62]. Park et al. [31] carried out an MEEFT cal- 
culation of S , Hep(0) with the use of da determined from (the same procedure as 
used in calculating S pp (0)), and obtained 

S Hep (0) = (8.6 ± 1.3) x 10- 20 keVb, (16) 

where the error spans the range of the A NN dependence for A NN = 500 - 800 MeV. 
Again, the MEEFT result agrees with that obtained in SNPA by Marcucci et al. [32] : 
5 , He P (0)sNPA = 9.6 x 10~ 20 keV-b. The above-mentioned large cancellation between 
the 1-body and 2-body contributions in this case amplifies the cutoff dependence 
of Snepify ) but the error quoted in eq.(16) is still small enough for the purpose of 
analyzing the existing Super-Kamiokande data [6]. 



8 The authors of this work use the term EFT* instead of MEEFT. 
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• Hen reaction 

The Hen reaction, 3 He +n — > 4 He +7, is very similar to Hep except, of course, 
the weak interaction is replaced by the EM interaction. Thus, just like the Gamow- 
Teller matrix element for Hep, the leading-order 1-body Ml matrix element for Hen 
is highly suppressed due to the pseudo-orthogonality of the initial and final nuclear 
wave functions. Moreover, there is strong cancellation between the hindered 1-body 
contribution and the 2-body exchange- current contribution, leading to a further 
suppression of the Hen reaction rate. Thus, like Hep, Hen is sensitive to higher 
order terms in chiral expansion. Meanwhile, the experimental value of <Ju e n is known 
rather accurately: cr^, = (54 ± 6) fib [63], or cr^ n = (55 ± 3) fib [64]. Therefore, 
Hen provides a good testing ground for the validity of MEEFT as used for Hep. A 
very elaborate MEEFT calculation of Hen was recently carried out by Lazauskas, 
Song and Park [65]. Up to N 3 LO considered by these authors, the Ml operators 
contain two LECs, which can be determined from the magnetic moments, /i( 3 H) and 
/i( 3 He). The initial and final nuclear wave functions are obtained from the rigorous 
Faddev-Yakubowski equations [66] for five sets of high-precision phenomenological 
potentials. Lazaukas et aVs calculation gives 

^Hcn = (38~58)/i6, (17) 

with high stability against changes in Ann in the range 500 ~ 900 MeV. The authors 
emphasized the importance of long-range nuclear correlation effects in obtaining an 
accurate estimate of o"H C n- To explain this point, it is useful to first discuss what 
Song et al. [67] pointed out in their MEEFT study of the Ml properties of the A=3 
nuclei. They noticed that the Ml matrix elements (MEs) in the A=3 nuclei cal- 
culated in MEEFT for various realistic potentials have strong correlation with the 
triton binding energies -B3 calculated with the corresponding realistic potentials. 
In fact, these quantities fall on a well-defined single curve. Meanwhile, since B 3 
governs the long distance contribution to the MEs, the model dependence (that is, 
variations in MEs calculated for nuclear potentials that give different values for i? 3 ) 
cannot be lifted by renormalizing the local (or short-range) operators. However, the 
use of the empirically found correlation curve between the MEs and B% allows us to 
drastically reduce scatters in the calculated MEs. This is achieved by introducing 
the constraint that one should accept only those values of the MEs which, along the 
correlation curve, have values of B3 consistent with its experimental value. This con- 
straint is found to essentially eliminate the model dependence in the MEs [67]. Now, 
Lazauskas et al. [65] noticed that a very similar feature exists also for Hen. Defining 
the quantity ( by ( = [q(a n H e 3 / /r H e 4 ]~ 2 ' 75 (*? = P n °ton momentum), Lazauskas et al. 
found that Ohcii and ( are so well correlated that the 0"H en -C plot forms a sharply 
defined straight line. Then, again, by accepting only those values of anon which are 
consistent with the empirical value of (, one can deduce the model dependence of 
UHen significantly. The estimate of ohcii given in eq.(17) was obtained through this 
method. 9 The fact that this theoretical value is in good agreement with the above 
quoted experimental value renders strong support for the basic soundness of the 

9 The main point here is not a linear relation between cthoii and £; after all, the latter is 
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MEEFT approach. We consider it noteworthy that MEEFT works well also for a 
highly suppressed process wherein higher-chiral order terms become important; see 
Ref. [68]. The examples of Hen and the MEs in the A=3 systems also demonstrate 
that it is profitable to avail ourselves of the merit of EFT (to control short-range 
physics) and the strength of SNPA (to control long-range nuclear correlations), as 
we do in MEEFT. 

3.5 ji - 3 He capture 

Ackerbauer et al. [69] carried out a high-precision measurement of the capture rate 
for /i~He 3 — > v^t. Achieving astonishing accuracy of 0.3 %, they reported the ex- 
perimental value r(/i 3 He) cxp = 1494(4) s _1 . A detailed MEEFT calculation of this 
capture rate was performed by Gazit [70]. Here again, the only unknown LEC, oIr, 
was fixed by fitting to for each given value of the cutoff parameter A NN (ranging 
500 ~ 800 MeV), and then r(/z 3 He) was calculated. The results are summarized 
as r(/i 3 He) = 1499(2)A(3)ATM(5)t(6)#c*, where the first error comes from the Ann 
dependence, the second error from uncertainties in extrapolating the nucleon weak 
form factors to finite momentum transfer, and in the choice of specific nuclear mod- 
els, the third error from experimental accuracy for Fg, and the fourth error from 
uncertainty in radiative corrections. The theoretical value (with total error ~ 1 %) 
agrees well with r(/i 3 He) cxp , providing a yet another case to support the validity of 
MEEFT; see also Gazit et al. [71], where a full EFT calculation for tritium /3-decay 
and the tritium binding energy is presented. 

In Ref. [59], quoted earlier in connection to S pp calculations, Marcucci et al. 
also carried out an elaborate calculation of r(/i 3 He) in both SNPA and MEEFT, 
and they found excellent agreement between the two approaches. Combining the 
results coming from these two methods, Marcucci et al. arrived at a highly accurate 
estimate T(/i 3 He) = 1484 ± 13 s -1 . Again, the theoretical value agrees very well 
(within quoted errors) with the above experimental value. 

3.6 Neutron-neutron scattering length a nn 

In the construction of a high-precision phenomenological NN potential, V N ^ nom , the 
experimental errors in the nn scattering length, a nn , are embedded in the overall x 2 - 
value, and it is not obvious to what extent the existing samples of V^ nom directly 
reflect the uncertainties in a nn . Meanwhile, low-energy electroweak transitions in 
lightest nuclei are in general rather sensitive to the precise value of a nn . For instance, 
in [4], uncertainty in the pp-fusion S -f actor, S pp (0), due to the existing uncertainty in 

a rather "esoteric" quantity, introduced just to "dramatize" the correlation between <7Hen and 
binding properties of the relevant nuclei. The message here is that, to get a reliable estimate of 
CHen, one must employ, even among so-called high-precision nuclear potentials, only those that can 
reproduce the binding properties of the A— 3 and 4 nuclei with needed accuracy. This criterion 
does not require the introduction of (. 
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a nn is estimated to be 0.5%, which is not devastatingly large but rather significant. 
Apart from this point, it is of general importance to determine, as accurately as 
possible, the empirical value of a nn , one of the fundamental parameters character- 
izing the NN interaction. Since extraction of a nn from the analysis of nd — > nnp 
reaction at present involves some inconsistencies [72, 73, 74], the standard value 
inn = — 18.63 ±0.27(exp) ±0.30(theor) fm [75] is based on the ir~d — > ^nn reaction. 
The theoretical framework responsible for this result uses phenomenological NN in- 
teractions, and intrinsic short-distance physics uncertainties in this framework limit 
the theoretical precision. Gardestig and Phillips [76, 77] investigated the possibility 
of using EFT to reduce theoretical uncertainties in the extraction of a nn . They have 
shown that short-range physics involved in the reaction can be parameterized using 
the same LEC, dn, as introduced in [31], and that, if for the sake of demonstration 
the currently accepted value of dn is used, the theoretical uncertainty in a nn can be 
reduced down to Aa„ n (theor) = 0.05 fm. This is quite an impressive reduction (by 
a factor of 6) of the existing theoretical uncertainty. 

3.7 Radiative corrections 

At the level of a few percent precision, radiative corrections (RCs) become an im- 
portant issue. For instance, in fip capture, the electroweak RC in the semileptonic 
reaction and the QED corrections to the muonic atom wave function are found to 
increase the capture rate by 2.8 ±0.4 % [78]. The best available estimates of RCs for 
single-nucleon processes have been given by Marciano and Sirlin [79] for neutron (3 
decay, and by Czarnecki et al. [78] for jj,p capture. In these works RCs of order a is 
decomposed into so-called outer and inner corrections [80]. The outer correction is 
a universal function of the lepton energy and is independent of the details of hadron 
physics, whereas the inner correction is affected by short-distance physics and the 
hadron structure. The inner corrections coming from 7- (weak boson) loop diagrams 
are divided into high-momentum and low-momentum parts. The former is calcu- 
lated in the current quark picture, while the latter is evaluated with the use of the 
phenomenological weak form factors of the nucleon. As regards RCs in multinucleon 
systems, all existing works (except those related to the + — > + Fermi transitions), 
concentrate on the single-nucleon RC diagrams because the genuinely two-body- 
type RCs are expected to be small. Kurylov et al. [81] carried out extensive studies 
of RCs in the A=2 systems based on the Marciano-Sirlin approach [79], and their 
results have been widely used in analyzing the A=2 weak-interaction observables 
(including those discussed earlier in this article). According to [81], with the use of 
Gf from /i-decay, and the "effective" g\ obtained from neutron /3-decay, radiative 
corrections specific to pp-fusion is estimated to be ~ 3-4 % (with 0.1 % uncertainty). 

Although the existing estimates of RCs based on the Marciano-Sirlin ap- 
proach [79] are believed to be reliable at the level of quoted accuracy, it is worth 
noting that HB^PT offers a model-independent framework to calculate RCs, insofar 
as all the relevant LECs are known. The first HB^PT evaluation of RCs for neutron 
f3 decay was carried out to order a by Ando et al. [82]. In this approach the short- 
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distance physics contributions can be condensed into two LECs, one related to the 
Fermi constant Gp and the other to the axial coupling constant g^- If one can deter- 
mine these LECs from empirical data, one would have a totally model-independent 
way to evaluate RCs to neutron (3 decay. A similar situation is expected to exist for 
any low-energy weak-interaction processes on a single nucleon. Even if one cannot 
determine these LECs at present, it seems of great interest to correlate RCs for var- 
ious weak processes on a single nucleon taking advantage of the fact that we know 
exactly what LECs are involved. 

4 Summary 

Despite the limited scope of topics covered, we hope we have succeeded in conveying 
the message that MEEFT is a highly reliable predictive framework for computing 
transition amplitudes for a large class of electroweak processes in light nuclei. We 
remark that, in each of the cases for which both SNPA and MEEFT calculations 
have been performed, it has been found that the results of the two methods agree 
very well. 

Although MEEFT has been registering and will continue to score significant 
successes, we should mention that great progress is being made in attempts to carry 
out full EFT (or "AB-EFT" ) calculations in which both the transition operators and 
nuclear wave functions are derived from HB^PT; see. e.g., [83, 84, 85]. We believe 
that, as increasingly higher order terms are included in AB-EFT calculations, their 
results will come increasingly closer to those of MEEFT, which are anchored on 
accurate experimental inputs. Although the developments in AB-EFT calculations 
are enormous and highly impressive, we will not discuss them here; we trust they 
will be covered in Ulf Meissner's contribution to this volume. 

As mentioned, Riska and Brown's ground-breaking work on the np — > 'yd 
reaction almost three decades ago was a precursor of explosive developments in 
high-precision calculations of electroweak processes in light nuclei. We hope Gerry 
is happy to see the enormous growth of the field in which he has played such an 
important role. 
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